Atomistic simulations of heat transport in real-scale silicon nanowire devices 
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Utilizing atomistic lattice dynamics and scattering theory, we study thermal transport in nanode- 
vices made of 10 nm thick silicon nanowires, from 10 to 100 nm long, sandwiched between two bulk 
reservoirs. We find that thermal transport in devices differs significantly from that of suspended 
extended nanowires, due to phonon scattering at the contact interfaces. We show that thermal 
conductance and the phonon transport regime can be tuned from ballistic to diffusive by varying 
the surface roughness of the nanowires and their length. In devices containing short crystalline wires 
phonon tunneling occurs and enhances the conductance beyond that of single contacts. 



Silicon nanowires (SiNW) can be now routinely synthe- 
sized via different processes, such as VLS or etching, with 
very good control on diameter, length, alignment, growth 
direction and surface properties^. Thus, they can be em- 
ployed in electronic applications as parts of solid-state 
memories^ or field-effect transistorsM. Even though bulk 
Si is a poor thermoelectric material, SiNW have also been 
proposed as promising building blocks for thermoelectric 
applications 5 . Relatively high thermoelectric figures of 
merit were measured at room temperature, mostly due 
to a large reduction of lattice thermal conductivity 6 -! 7 -. 
There is indeed experimental evidence 8 - - — , supported by 
theoretical studies^— that the thermal conductivity of 
silicon at room temperature can be reduced by as much 
as two orders of magnitude by nanostructuring and sur- 
face treatmen t 16 ' 17 . 

For any technological purpose it is essential to control 
heat transmission and dissipation in the device, carefully 
resolving the contribution from contact interfaces and/or 
supporting substrate s 18 ' 19 . Single crystalline SiNW with 
monolithic lead contacts were recently fabricated using 
electron beam lithography^. Their measured thermal 
conductivity (20 Wm^'K" 1 ) resulted significantly higher 
than that of previous experiments on SiNW of similar 
diameter— (2 Wm _1 K _1 ), indicating that contacts play a 
major role in determining the overall thermal properties. 

So far, atomistic modeling has been limited to ex- 
tended periodic systems, whereas simulations of devices 
are commonly performed by mesoscale models^—. Here 
we utilize atomistic lattice dynamics to investigate ther- 
mal transport within SiNW devices made of 10 nm thick 
and up to 100 nm long wires connected coherently to two 
semi-infinite bulk reservoirs. The device thermal conduc- 
tance is computed for wires of different lengths, with ei- 
ther crystalline or amorphized surfaces. We discuss how 
length, contacts and surface disorder determine different 
transport regimes. In particular, our calculations show 
that thermal transport in sub-micrometer devices is in- 
herently different from that of extended wires, as it is 
controlled in large parts by contacts and finite size of the 
devices. 



The wires are grown in the (100) crystallographic direc- 
tion, have a nearly circular cross section and form ideal 
monolithic contacts with the reservoirs, preserving stack- 
ing and crystallographic orientation. Also, the surface of 
the crystalline SiNW is reconstructed in order to mini- 
mize the number of dangling bonds^. Atomic positions 
of the wires are relaxed in a periodic cell at zero tem- 
perature and the cell length optimized so to obtain an 
unstrained system. The reservoirs are considered peri- 
odic in the xy plane and semi-infinite in the z plane. For 
both relaxation and heat transport calculations, the in- 
teractions between silicon atoms are modeled by means 
of the short-range empirical forccficld by Tcrsoff 25 . 

We compute the thermal conductance of SiNW de- 



30nm 




FIG. 1. Representation of the core-shell silicon nanowire de- 
vices wires of different lengths, for which we have computed 
thermal conductance and heat transport properties. The sil- 
icon nanowires are connected coherently with semi-infinite 
bulk silicon thermal reservoirs. The core-shell nanowires 
shown here have a crystalline core of ~ 8 nm in diameter 
and a ~ 1 nm thick amorphized shell. Surface amorphization 
is obtained by heating the system slightly beyond its melting 
temperature (2500 K for the empirical model of silicon em- 
ployed here) for 300 ps, quenching to room temperature in 
700 ps, and eventually optimizing the structure by steepest 
descent. 
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vices by harmonic lattice dynamics^—, employing the 
scattering matrix approach 31 . This method permits the 
calculation of quantum thermal transport in open sys- 
tems, yielding transmission spectra, thermal conduc- 
tance and real-space representation of heat flux and non- 
equilibrium energy density. The latter is obtained by 
computing the local density of states (LDOS) and allows 
us to probe different transport regimes. In spite of its 
assumptions and approximations, the most important of 
which is considering only elastic scattering, it has been re- 
cently demonstrated that this approach can consistently 
reproduce the Kapitza resistance of a silicon/germanium 
interface computed by molecular dynamics up to fairly 
large temperatures (T ~ 500 K) 32 . At higher temper- 
atures anharmonic effects should be explicitly consid- 
ered^ 3 .. 

a. Crystalline wires We first investigate the effect 
of varying the length of crystalline SiNW between two 
bulk reservoirs. Wire lengths of 10, 20, 30, 50 and 100 
nm are considered: the largest system is made of 6 • 10 5 
atoms. The transmission spectra are shown in Fig. [2^, 
and compared to the transmission spectrum of the in- 
terface between semi-infinite bulk and the semi-infinite 
SiNW2i. 

We find that the transmission curves overlap for fre- 
quencies up to 4 THz regardless of the wire length 
(Fig. indicating that phonon transport at low fre- 
quency is ballistic. In this frequency range, there is no 
attenuation within the wire itself but only scattering at 
the bulk-wire interface. The device behaves as a complex 
Fabry-Perot cavity where the large number of channels 
average out the interference patterns. As a consequence, 
the transmission is independent on the length of the de- 
vice, and so is the conductance, er(T), in the low temper- 
ature range. The ballistic regime below 4 THz was also 
probed computing the LDOS (Fig. [3^), which indeed dis- 
plays no attenuation within the wire. The LDOS decays 
sharply at the contacts and oscillates around a constant 
value. Those oscillations are reminiscent of the interfer- 
ence patterns formed within the device acting as a cavity. 

Surprisingly, the transmission of the devices becomes 
larger than that of the semi-infinite wire 3 -^ for acoustic 
frequencies larger than 4 THz. We observe that the to- 
tal transmission may even exceed the number of available 
channels in the wire (Fig. 4), indicating that a number of 
bulk states tunnels through the device. As a result, the 
total thermal conductance (Fig. [2p) turns out larger than 
that of the single contact for wires shorter than 30 nm. 
The contribution from the tunneling states to a decreases 
when the length of the wires is increased, and vanishes 
for 100 nm and longer devices. Plotting the local energy 
density allows direct visualization of the tunneling trans- 
port regime (Fig. [3}}). In contrast with the low frequency 
ballistic modes, the LDOS does not drop abruptly at the 
interface, but it penetrates the device region and decays 
slowly. A natural consequence of the combination of bal- 
listic transport and phonon tunneling is that the thermal 
resistance of the two contact interfaces is not additive 
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FIG. 2. Transmission spectra (a) and integrated ther- 
mal conductance (b) of devices containing crystalline silicon 
nanowires with 10 nm diameter and length varying between 
10 and 100 nm (dashed lines), compared to the transmis- 
sion spectrum and thermal conductance (a(T)) of a single 
bulk/nanowire contact (solid line). The inset in panel (b) 
zooms in the low temperature region of cr(T). 



even for wires as long as 100 nm, contrary to what one 
would get for diffusive transport and contacts in series. 

We also compare our computed thermal conductance 
with experimental measurements at very low temper- 
ature^, though they refer to much thicker nanowires 
(~ 100 nm). In the temperature range between 5 and 
20 K, cr(T) is proportional to T 3 , (inset of Fig. [2p), in- 
dicating that dimensionality reduction has no significant 
effect on the low temperature trend of a. The T 3 trend is 
the same as found for single contacts 3 ^, suggesting that 
at low temperature heat transport in these devices is de- 
termined by scattering at the contact interfaces. At vari- 
ance with Ref£ we do not observe a temperature range 
where <r(T) oc T 2 when SiNW are crystalline. This dis- 
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FIG. 4. Difference A between the transmission spectra of de- 
vices containing crystalline silicon nanowires and the trans- 
mission spectrum of an infinitely long silicon nanowire of the 
same type and diameter (i.e. its number of conduction chan- 
nels). We consider devices made of wires of various length 
between 10 an 100 nm and diameter of 10 nm. For devices 
containing short wires (< 50) nm the difference is positive 
in a wide range of frequencies larger than 4 THz, indicating 
that in this range of frequency more modes are transmitted 
across the devices than the available conduction channels in 
the wire. Some bulk modes can tunnel through the device. 
The physical principle for phonon tunneling is analogous to 
the well known model of quantum tunneling across a finite 
barrier. The difference between the transmission spectrum of 
the contact interface and the number of channels of 10 nm 
wire is also reported for sake of comparison. 



FIG. 3. Local density of states of crystalline wires (a-b) and 
core/shell wires (c). Each representation is associated to left 
hand incoming phonons contributions for the specified fre- 
quency. 

crepancy is most likely due to the perfect crystallinity of 
our models, and different trends are observed when the 
surface of the SiNW is amorphized. 

b. Core-shell wires We now discuss thermal trans- 
port in devices with crystalline core / amorphous shell sil- 
icon nanowires (CS-SiNW). To obtain CS-SiNW, crys- 
talline SiNW are heated beyond the melting temperature 
of the bulk for a short time (300 ps) , quenched to room 
temperature in 700 ps and optimized using the steepest 
descent algorithm. This procedure produces a ~ 1 nm 
thick amorphous layer. Four models have been designed 
for this purpose: two with 30 nm long wires, one with 
60 nm and one with 90 nm long wires. The surface dis- 
order in the first 30 nm sample is limited to about two 
surface layers, whereas the amorphous shell of the re- 
maining 30 nm, 60 nm and 90 nm samples is about 1 nm 
thick (Fig. 1). The transmission spectra and the thermal 
conductance of these devices are shown in Fig. [5] With 
respect to crystalline devices with wires of similar length, 



we observe that shell amorphization leads to a reduction 
of the device thermal conductance by a factor between 3 
and 6 at room temperature. The larger reduction occur- 
ring in the longer device. For sub-micrometer wire de- 
vices with diameter above 10 nm, surface amorphization 
does not have an effect on thermal transport as dramatic 
as that observed in thinner SiNW—. This indicates that 
when large thermal resistance is experimentally observed 
for SiNW thicker than 10 nm, surface disorder has to 
be supplemented by other sources of phonon scattering, 
such as extended defects or poor contact interfaces. 

The inset in Fig. [5Jd shows that at low temperature 
a(T) oc T 2 ' 2 . This indicates that elastic scattering at the 
crystalline/amorphous interface can be the cause for the 
T 2 trends of thermal conductance observed experimen- 
tally in the temperature range between 1 and 10 K— . Di- 
mensionality reduction plays the indirect role of enhanc- 
ing surface effects, due to the higher surface-to-volume 
ration in nanostructures. At variance with experiments, 
the T 2 range for cr(T) extends to 50 K: this discrepancy 
may be due either to size differences between models and 
experiments (experimental SiNW have much larger cross 
section), or to neglecting of anharmonic scattering pro- 
cesses in our lattice dynamics approach. 

The calculation of LDOS in the CS-SiNW devices 
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FIG. 5. Transmission spectra (a) and integrated thermal 
conductance (b) of crystalline core-amorphous shell silicon 
nanowires with 10 nm diameter and lengths varying between 
30 and 90 nm. The inset in panel (b) zooms in the low tem- 
perature region of <j(T). 



shows that phonons with frequency up to 1 THz are 
transmitted almost ballistically along the wire. The 
transmission coefficient is not sensitive to either the 
wire length or the thickness of the amorphous layer, 
and is comparable to the one of the crystalline devices. 
(Fig. [5^). At higher frequencies, the phonons undergo 
multiple scattering from the rough crystalline- amorphous 
interface and the LDOS decays linearly along the wire, 
as shown in Fig. [3p. The transmission coefficient drops 
accordingly, becoming significantly lower than the one 
of the crystalline device. This behavior is characteristic 
of the diffusive transport regime. We also observe that 
the LDOS is larger in the amorphized shell of the wires, 
which seems to indicate that the shell is warmer than the 
core in non-equilibrium conditions. For phonon frequen- 
cies up to 2 THz, the LDOS and heat flux retain similar- 
features as for crystalline wires. At higher frequencies, 
the heat flux is concentrated within the core region of the 
wire and presents irregularities, confirming the effect of 
the shell amorphization on the wire conduction proper- 
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FIG. 6. Ratio between the resistance of a single crystalline 
contact and the total resistance of the devices with core-shell 
nanowires 30, 60 and 90 nm long. 

Assuming that the contact resistance does not vary sig- 
nificantly from crystalline to core-shell wires, we use the 
values obtained in a previous workai, and compute the 
ratio between a single contact resistance and the total re- 
sistances of the core-shell SiNW devices (Fig. [S]) - The rel- 
ative importance of contact thermal resistance depends 
on the temperature and on the length of the wires. Con- 
tact resistance prevails at low temperatures, reaching a 
maximum over 60% of the total device resistance below 
10 K, and becomes less important at higher tempera- 
tures. In the low temperature regime only low frequency 
modes are populated, which conduct ballistically through 
the wires: contacts are the only source of scattering. For 
T larger than 100 K the curves in Fig. [5] reach a plateau, 
which value depends on the length of the wires. The rel- 
ative effect of contact resistance at T > 100 K decreases 
from ~ 35% for the 30 nm wire device to ~ 15% for the 
90 nm wire device. For longer wires the relative weight 
of contact resistance would decrease further, as the wire 
resistance would increase linearly with the length, yet 
our estimate is a lower boundary, since we consider ideal 
monolithic contacts. For instance, if we assume that we 
can treat the 90 nm wire and the contacts as thermal re- 
sistors in series, and evaluate each contribution as addi- 
tive, then the contacts contribute about 30% of the total 
thermal resistance of the device at room temperature. 

In summary, our calculations prove that: i) Ther- 
mal transport is ballistic in devices containing crystalline 
nanowires, and scattering of acoustic phonons occurs 
mostly at the contact interfaces. However, when SiNW 
are shorter than 30 nm the thermal conductance may be 
larger than the contact conductance, due to tunneling 
of bulk vibrational modes through the SiNW. ii) Amor- 
phization of a surface layer can reduce the thermal con- 
ductance of these type of devices by at most six times at 
room temperature, as it is ineffective in scattering low fre- 
quency acoustic modes. The introduction of long-range 
disorder or nano-structuring is probably essential to fur- 
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ther decrease the low frequencies contribution to thermal 
conductance, iii) The contribution of contact interface 
scattering to the thermal resistance is significant even 
in those systems where there is strong surface scattering 
and heat transport is mostly diffusive, iv) In CS-SiNW 
devices the thermal conductance at low temperature is 
proportional to T 2 , in remarkably good agreement with 
experiments^. We note that the surface scattering term 



is sufficient to account for this unusual trend of cr(T). 
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